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Abstract 



The evolution of inhomogeneities in a spherical collapse model is studied 
by expanding the Einstein equation in powers of inverse radial parameter. 
In the linear regime, the density contrast is obtained for flat, closed and 
open universes. In addition to the usual modes, an infinite number of new 
growing modes are contained in the solutions for pressureless open and closed 
universes. In the nonlinear regime, we obtain the leading growing modes in 
closed forms for a flat universe and also, in the limits of small and large times, 
for an open universe. 



I. INTRODUCTION 

The Newtonian and the relativistic theories of small perturbations form the theoretical 
basis for our present understanding of the formation of structures in the universe [0-0. In 
the Newtonian regime, which is appropriate for pressureless universes, the hydrodynamic 
equations, incorporating the expansion of the universe, are perturbed by small fluctuations. 
Subsequently, a second-order differential equation is obtained whose solutions give the de- 
caying and the growing modes of the density contrast. This quantity measures the evolution 
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rate of the perturbation with respect to the background density and thus determines whether 
structures can evolve from small initial perturbations or not. 

In the relativistic regime, one expands the Einstein metric g^ v around the Friedmann- 
Robert son- Walker metric by a local perturbation, h^ v , which arises naturally in an imperfect 
fluid model of the universe. Consequently, the evolution of the perturbations is formulated 
by means of the energy conservation equation, which can be written in terms of the time 
derivatives only, and a combination of the linearized Einstein equations (see e.g. 0, chap- 
ter 15). Space derivatives are eliminated in this procedure and a second-order ordinary 
differential equation in time is obtained for the density contrast. 

The procedure in the relativistic regime can be summarized as follows. Fixing the gauge 
such that the perturbation vanishes in the (iO) and (00) directions, one obtains, for the 
space-space component of the Einstein equation, 

\7 2 h - d 2 hik _ d 2 hjk d 2 hkk _ ^>2^ 
13 dxWx k dx i dx k dx % dxi 13 

+RR (hij - Stjhkk) + 2R 2 {-2hij + Sijhkk) = (1) 
— 8nG (p — p) R 2 hij — 8irGR Sij (5p — dp) + ( imperf. fl. correc.) , 

where the last term corresponds to the imperfect fluid corrections which are given in terms 
of a velocity field. The differential equation contains space derivatives of the metric and, 
therefore, is extremely difficult to solve. In general, matter distribution constrained to a 
small region has large space derivative as compared to the metric itself. This is clear for 
a well-localized Gaussian distribution, p(x) ~ e -K x ~ x o) 2 ^ f or j ar g e \^ which cannot be ex- 
panded in inverse powers of the distance rajgn ■ The off-diagonal components of the Einstein 
equations are not amenable either, once again, due to the space derivatives. However, the 
time-time component contains no space derivative and reads 

hk ~ 2^h kk + 2 ((^j - | j h kk = -8nG {8p + 3^) R 2 . (2) 

Additional informations can also be obtained from the energy-momentum conservation equa- 
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tionQ, 

$P + 3^ (6p - dp) = (P-P)^ f + ( im P erf - fl - correc.) (3) 

which, together with a few general hypothesis, leads to the simple relation 5(t) = 8p/ {p+p) = 
—h} t k/2R 2 . This relation, when used in (|2|), gives rise to a second-order differential equation 
whose solutions give the growing and the decaying modes of the density contrast S(t). 

In the nonlinear regime, evolution of mass density fluctuations has been studied in the 
matter-dominated Einstein-de-Sitter universe by means of the power spectrum, which is 
obtained from the Fourier transform of the density contrast |4j] . 

Although the standard procedure is formulated in such a way as to enable us to describe 
the evolution of a general kind of perturbation, the hypothesis of an imperfect fluid with 
a velocity field which is at the root of this scheme cannot be applied to certain kinds of 
fluctuations. In the standard scheme, the fluctuation seeds, which are believed to give 
origin to the galaxies, are described by Poisson-like perturbations. These are well-localized 
functions with large space derivatives. 

Contrary to the standard scheme, we study soft perturbations, which are almost homo- 
geneous at large scales. These perturbations are especially interesting for the purpose of 
studying inhomogeneities at very large scales [fH|. They are also important for studying 
totally inhomogeneous universes UTHlOU . Although perturbative expansion around a comov- 
ing metric is a subtle procedure, smooth perturbations at large scales can be incorporated 
into perturbative schemes. 

We consider the homogeneous Friedmann- Robert son- Walker metric as a background met- 
ric and expand the perturbation in powers of the inverse radial parameter, i. e. in powers of 
1/r. The advantage of this procedure, which is valid for very large scales, is that we can con- 
trol the space-dependence of the perturbation at each order, and, as we shall see, the Einstein 
equations will provide us with more informations on the evolution of the perturbations. 



x It is worth commenting that the inclusion of higher derivatives may introduce zero modes. 
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In this approach, equation (P simplifies and all the space derivatives disappear in the 
first approximation (the zero-th order terms, being homogeneous, are zero modes of space 
derivatives). In particular, even though (Q) is still valid, we no longer need higher-order time 
derivatives to determine the density contrast. 

Our aim is to use the aforementioned expansion in 1/r for spherically symmetric config- 
urations, i.e., for the S- waves, and verify the presence of the growing modes for the density 
contrast. 

We parameterize the metric as 

= dt 2 _ R ^ r)2 + r w) (4) 

where the scale factor Rp(t,r) is a function of the radial parameter as well as time. Using 
this metric in the Einstein equation we obtain 

.ftp y rip J JXp Tip 

-^(l + 2 W )^ + (l-.)^ + 2^)(l-^ 2 )=0. (5) 

where "dot" and "prime" stand for time and space derivatives, k = 0,1, —1 for flat, closed 
and open universes respectively and the parameter u is the pressure-density proportionality 
constant, i.e. P = up. The density is given by the time-time component of the Einstein 
equation 

3 iw) 2 + + §C-2^4) (1- + J* + 2kr^ = 8,G Pp (r, t). (6) 

\ n P / \ r -Kp Tip ftp J Tip Tip 

To solve the above equations we expand the scale factor and the density around the 
Friedmann background by small fluctuations in powers of inverse radial parameter. That is, 

oo 

R p (t,r) =R+Y / $Rn(t)r~ n , (7) 

n=0 
oo 

p P (t,r)=p + Y / Sp n (t)r- n , (8) 

n=0 

where R and p are the Friedmann scale factor and density respectively. It is worth com- 
menting that these solutions are not just gauge modes of the metric but represent genuine 



perturbations |. It is easy to see that a solution of 5g^ u = + £ VjlM where 5g^ u is the 
metric perturbation, such that 5go r = and 5g rr ~ t a r~ n , cannot be obtained. 

In the following sections, we solve the perturbed Einstein equations for different values 
of k and u and find the density contrast which is defined as: 

*,0 = ^ (9) 
where 5p is given by the summation term in equation (Rp. 

II. MATTER-DOMINATED OPEN UNIVERSE 

The scale factor in the Einstein equation fl5|) is expanded as in (0) and the coefficients at 
the r~ N order are collected. For u = and k — — 1, these coefficients satisfy the recursive 
relation 

^ R n 5 R m 5 Ri5 FiN-n-m-l + &Rn&Rm&Rl&RN-n-m-l) + 

n,m,Z 

+ ^ n(4-iV + ^^^-3-^ = 0, (10) 

n 

where 5i? n represents R + 5Rq for n = and the indices n, m and / run from zero to infinity. 
The above equation can be solved recursively for each N. In general, using the parametric 
expressions R = cosh^ — 1 and t = sinh'0 — ip, we obtain an inhomogeneous second- 
order differential equation, whose non- homogeneous terms are given in terms of 5Rj <n (ij)). 
Although, it is not possible to solve such an equation in the general case, the homogeneous 
part of equation (|10"D , 

d 2 SR n ( 1 



# 2 ' \ U 2sinh 2 f, 



SR n = , (11) 



can be solved recursively for each n. The corresponding solutions are 



2 For a comprehensive study of the gauge freedom arising in the context of density perturbations, 
we refer the reader to reference JL 1[J . 
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(5R n ) + ~ —= — — -— — (cosh.(if)/2) sin \fnif> — 2y/nsmh(if)/2) cos y/n 
ynsinh(^/2j v 

(5R n )„ ~ —r, /o\ (cosh(V>/2) cos + 2Vrisinh(V>/2)sin > ( 12 ) 

smh(ip/2) v ' 

Using this result in (^), we obtain the density and subsequently the density contrast by 
using expression (0). This is given by the growing mode 



n 



and the decaying mode 



3 sinh ip 2n sinh if) 

l-cosh?/>) 2 1-cosh^ 



, , / \ 3 sinh if) 2n sinh if) , / x. o -|- cosn 

L w ~cos vnw rr H ; — - — vnsm u/nw — - . (14) 

VV ^[(l-cosh^) 2 1-cosh^J V ^l-cosh^ v ; 

The n = (ie. the r-independent perturbation) expressions reproduce the standard growing 
mode 

+ (l-cosh^) 1-cosh^ ' v ; 

and the standard decaying mode 

«T-W ~ 77~~~T~7T2 (16) 
(1 — cosh-?/;) 



+cos (^)r^h^ ' (13) 



5 + cosh if) 



for the density contrast JT|Q|. The solutions ( 13 ) and flT4] ) show that in addition to the 



standard growing mode, infinitely many growing modes exist at all higher orders in the 
perturbation expansion in the linear regime. 

All the above results can be used for the closed universe by making the substitution 
if) = -i9. 

Non-linear growing modes would be expected if one could solve the inhomogeneous dif- 
ferential equation flTOl). In Section V, we obtain these modes in the two extreme limits of 
small and large times. In the small-time limit, open and close universes behave as a flat 
universe. In the asymptotic, i.e. large time, limit the scale factor varies linearly with time 
and the inhomogeneous equation (|1(]) can once again be solved at higher orders. We obtain 
a new growing mode in this limit. 
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III. MATTER-DOMINATED FLAT UNIVERSE 



For the flat universe, the series expansion (0) of the Einstein equation (§) leads to the 
recursive equation 

J2 (25R n 5R m 5Ri5R p + 5'R n 5'R m 5Ri5R p ) r ^ n+rn+l+p) 

n,m,l,p 

+ 0(2 - m)5R n 5R m ) r ^ n+m+ ^ = (17) 



where the indices n, m, I and p run from zero to infinity. Unlike, the differential equation 
(|lOl) for the open universe which could not be solved without discarding the inhomogeneous 
terms, the above inhomogeneous equation for the flat universe can be solved recursively to 
find the perturbation coefficients R[n](t) 0. 
The first few coefficients are 
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5R = c t 3 , (18) 
5R 1 = Cl t^ + c 2 t~^ , (19) 
SR 2 = c x ti + c 3 H - j4t~^ , (20) 

9 4 2 9 1 

= ~ io Clt3 + c ^ 3 + 2° 2t3 

1 /l 9 1 \ 4 1^7 . 

+c 5 t~ + y-CK? - -c 2 c 3 J rs + , (21) 

/9 2 9\ 4 2 

5^4 = ( jc( - -J ts + c 4 t» + 

99 \ i i 27 2 2 

2 

8 2 

" „ 10 



-yCl + 9c 3 J ts + C 6 t 3 + —eft 

+c M r* + c a H + C b rf , (22) 



243 2 / 27 81 o 36 2 \ 4 243 2 

5i? 5 = c x t + c 4 c? H cf t3 c 2 t + c 7 t3 + 

280 V 10 20 1 5 V 8 

153 333 27 513 2 \ 1 _i 

-C1C2 T7T c i c 3 + + —7-0^2 t 3 + c 8 £ 3 + 



5 10 2 10 

693 9 45 \ 2 ^ 4 _ 7 10 . . 

Cl c 2 + — c 2 c 3 rs + C c rz + C d t"3 + Cft" , (23) 



40 " 4 



3 The results have been reproduced by a Maple V-5, program 15] 
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where C"s are functions of either preceding c's or are new constants themselves. In general, 
the leading growing modes can be written in a closed form as 

oo j.2(n+2)/3 

^ + (M) = E4„ + 3-^T • ( 24 ) 

n=0 

From the time-time component of the Einstein equation (|6|), we recover the mass distri- 
bution and subsequently the density contrast 5(t,r) = J2^L ^n(t)r~ n , whose first few terms, 
up to an overall constant factor of 3/4, are 

tfi = -A^r 1 , (25) 
5 2 = 4 (cica - c 3 ) t" 1 + 9c 2 t~ 2 , (26) 

12 2 _1 

d 3 = — -cits - 6c 2 t 3 
5 

+4 (cic 3 + cic 2 + c 2 c 2 - c 5 ) + 18 (c 2 c 3 - cic 2 ) r 2 - 18c^" 3 , (27) 

/ 14 47 9 \ 2 /184 \ i o 4 

5 4 = ( - —ci + —c\\ t 3 + ( — cic 2 - 16c 3 J fi + 17c 2 r 3 

4 (c 2 c 4 + cic 5 + c 3 c 2 - c 6 + cic 3 - 2c 2 c 2 - c 2 c 3 ) + 

9 (-4cic 2 c 3 + 2c 2 c 5 - 2 Cl c 2 2 + 3c 2 c 2 2 + cfj t~ 2 + 54 ( Cl c 3 2 - c^c 3 ) t~ 3 + — c 2 t" 4 , (28) 
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162 4 / 81 , 

(5 5 = — cits + I -54c 2 - yC 2 ) ts 



+ -12c 4 - 4c 3 + 16c 2 + —c\ - —c\ - — c 4 ta + C gg t-3 + C kk t—i + C h t~ 
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132 „ 96 o 36 \ 2 „ _i 

— c? c? 

5 1 5 1 5 

ar 1 + c,t- 2 + cy~ 3 + ar A + c m r 5 . (29) 

In general, the leading growing modes can be represented in the closed form 

00 i2(n+l)/3 
fl+ = E4n+3 r2n+ 3 • ( 30 ) 

These results show that up to the third order in the perturbation expansion, the density 
contrast does not grow. This is due to the fact, mentioned in the introduction, that the 
perturbation coefficients at these orders are not sensitive to the space-dependence of the 
perturbed scale factor R p (t,r). However, the third-order perturbation term contains the 
usual growing mode in the linear regime Q, namely S ~ t 2 ^ 3 . (This term emerges because of 



It is worth making the comparison with the open universe where the standard zero modes 
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the t 4//3 growing mode in (21).) The nonlinear modes are given by n)0 terms in expression 
(|30|) and fully agree with the results of reference |4j]. 

Once again, we see that in addition to the usual growing modes arising in the linear 
regime, many more growing modes appear at higher orders in the perturbation series. The 
higher growing modes can have interesting consequences for structure formation since differ- 
ent perturbative regimes can correspond to different orders in the perturbation series. For 
very large structures, such as superclusters, the dominating terms are the first few, which 
vary with a smaller powers of time and so are expected to appear later. On the other hand, 
small fluctuations are described by higher-order terms in the perturbative series and thus can 
form earlier. This simple argument points toward a bottom-up scenario for the formation 
of structures and thus a universe dominated by Cold dark matter 



IV. RADIATION-DOMINATED FLAT UNIVERSES 



In this case, we solve, using a Maple V-5 program [15], a recursive equation similar to 
([HI) for the coefficients of the scale factor. As in the previous cases, the leading growing 
modes can be written in the following closed form: 



oo / rz\ 2 "+ 3 

5R + (r,t) = J2 c 2n + s — • (31) 

n=0 V r / 

Substituting this solution in the time-time component of the Einstein equation (Q) and 
subsequently using (fl we obtain the expression 

oo 

^(M) = E^ + 3-^ (32) 

n=0 ' 

for the leading growing modes of the density contrast. 

Once again, in addition to the standard growing mode, 8 ~ t, which arises at the third 
order of the series, we obtain an infinite number of higher-order growing modes. 



already appear at the zero-th order. 
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V. NONLINEAR MODES IN THE OPEN UNIVERSE 



In Section II, we obtained the density contrast for a pressureless open universe in the 
linear regime. We then commented that the nonlinear modes can be obtained in the asymp- 
totic limit and in the limit of small times. In this section, we study an open universe in 
these two extreme limits. Our analyses are not restricted to a dust universe. 

It is a non-trivial task to define the asymptotic domain in the open Friedman universe. 
Indeed, in an asymptotic limit compatible with the Einstein equation and with any linear 
equation of state, the scale factor increases linearly with time. This limit can be used, 
formally, to describe the asymptotic behaviour of the scale factor in the open universe, 
keeping in mind that, physically, the notion of large times is only meaningful when used 
relative to the Hubble time. On the other hand, observations are performed for small times, 
namely for the primordial universe, when open and closed universes behave as a flat universe. 

In this section, we discuss these two limits of our solutions, i.e. R ~ t, and R ~ t 2 / 3 . 
In the former case, the suppression of the growing modes implies that the inhomogeneities 
cease to grow after a sufficiently long time. Indeed, solving the series equation (|T0| ) in the 
asymptotic limit gives |I5| 



5R { 



o 



t 



(33) 



5R 



i 



Ci cos In t + C2 sin In t 



(34) 



SR 2 




68 



+few similar terms. 



(35) 



(36) 



Substituting the above solutions in the time-time component of the Einstein equation (^|), 



we recover the mass distribution and subsequently the following density contrast: 
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51 = 2ci cos In t + 2c2 sin In t — 6ci sin In t + 6c 2 cos In t , (37) 

5 2 = -Q^c 3 sin y/2l nt + 6\/2c 4 cos lnt + few similar terms. (38) 

The terms cos lnt and sin lnt do not represent any growing modes and only a decaying 
mode of 1/t is indicated by the above solutions. In the radiation-dominated era, however, 
we do find a growing mode. In this era, the first few corrections to the scale factor are 

5R = t , (39) 

1 \fl 1 v^39 
SRi = —7=Ci sin — In t H — =c 2 cos In t , (40) 

\/t 2 \/t 6 

1 /t 1 /7 

5^2 = ^c 3 cos^-lnt + ^c 4 sin^-lnt + C (r 3 ) , (41) 

2 yt 2 v 7 

and consequently to the density contrast are 



v 39 v 39 

^ = _[ Cl + v^csjv^sint^— lnt] + [-c 2 + civ / tv / 39]cos[^— lnt] , (42) 

6 6 

8 2 = — jj/s(0 sin[-^lnt] - 3c 3 Vt cos[-^ hit] - 3c 4 v / tsin[^ lnt] + j^U(t) cos[^lnt] 

. .^9, , r-s/39, , 7+131 2 ,^9, . 7+131 2 . 2r V39, . 

— 3ci Co sm lnt cos lnt N c 9 cos lnt N c 9 sm lnt 

6 JL 6 J 4t 2 1 6 J 4t 2 L 6 J 



5 /7TT-1, a 2x r v39 , . . r v39, . 
+ -v / 39-(c 2 -c 2 )cosP— lnt smp- lnt 
lib o 

1 /t 1 /t 

+^Vtf 4 (t) sin[— lnt] + -Vt/aW cos[— lnt] 

-3v^^c 3 sin^lnt-c 4 cos[^lnt]^ (43) 

where / 3 and / 4 are functions depending on integrals of the functions already appearing 
above. That is, 

/a = /25t- 5 / 2 c 2 cos(^lnt^ , (44) 

h = J 25t~^cj sin In t) - sin In tj cos 2 In tj + • • • . (45) 

We find the growing mode 5 ~ y/t already at the first order. Thus, asymptotically, the 
inhomogeneities can grow in an open universe during radiation-dominated era but stop 
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growing during the matter-dominated era. We interpret these results as a freezing of the 
present mass distribution []. 

At small times, the parametric solutions for the Friedmann scale factor in the matter- 
dominated open universe approach the solution for the flat universe, i.e. R ~ t 2 ^ 3 . However, 
in the above perturbation scheme we cannot simply abandon the constant k and approach 
the result from the point of view of the Einstein equations for a flat universe, since there are 
important contributions from the space-dependent derivatives through k (e.g. in equation 
(H)). Using again a Maple V-5 program [ITS] we find 



SR 


= cot% 


8R 1 


= citi + 


5R 2 


= c 2 P + 


SR 3 


= C 3 t3 + 



(46) 
(47) 
(48) 
(49) 

Substituting these in the time-time component of the Einstein equation @, we recover the 
mass distribution and subsequently the following perturbative contributions to the density 
contrast: 

s = c n t% , (50) 

51 = Cot" 1 + C p *f + • • • , (51) 

5 2 = C q r l + C s r 2 + C tt t! + . . . , (52) 

5 3 = c v r l + c w r 2 + c x r 3 + c v t* + c z t~* + c aa t 2 + ■■■ , (53) 

where C's are the constant coefficients. Two points are worth mentioning. In the open 
universe, growing modes appear earlier due to the r factors in the time-time component of 
the Einstein equation (Bl). In addition to the usual growing mode, namely t§, there exist new 



5 Interestingly, the seizure of the inhomogeneities to grow after a long time is compatible with 
the numerical results of reference [14| for a self-similar Newtonian universe. 
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growing modes which appear at lower orders as compared to the flat universe (see equation 

Moreover, later in the expansion of the scale factor the appearance of new growing 
modes due to the inhomogeneous part of the Einstein equation as written in (|T0| ) for the 
matter-dominated era proceeds as in the flat universe, but now with enhanced consequences 
for the density contrast (namely they appear at lower orders). This profile could not be 
studied directly in the parametric solution for the matter-dominated open universe due to 
the complications arising from the inhomogeneous part of the differential equation (|T0|) . 

In addition to the usual growing modes, the density contrast has an infinite number 
of higher-order growing modes. These are not expected in the usual linear perturbative 
schemes [TT| — |3|,p73 j . The higher growing modes were overlooked in the exact solutions for 



the open universe (|T2|) because only the homogeneous part of the differential equation was 
considered, while the new growing modes seem to be due to a back reaction of the system 
to the perturbation itself. 

The modes are divided in families, once they appear at some order of the perturbation, 
they continue appearing at higher orders together with their descendants which are one 
order smaller in time. 



VI. SMALL-DISTANCE PERTURBATION 

In the preceding sections, we have considered the perturbative expansion appropriate for 
large distances. We now analysis the perturbative expansion in positive powers of r which 
is more appropriate for small distances and thus more relevant to the formation of galaxies 
and other smaller structures. We start with the expansion 

oo 

B T (t,r) = R(t) + 1 £8Rn(t)r n , (54) 

n=0 

which is substituted into the Einstein equation ([5]). At order r -2 , the equation is trivially 
satisfied and the order r _1 terms imply that 5Ri vanishes. The next term, however, involves 
both R and <5i?2- If we insist on using the homogeneous result for R then the whole series 
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disappears. That is to say 8R n = for n 7^ 0. This problem can be understood once one 
realizes that there exists the simple separable solution 

** r > = TW ' (55) 

which implies that SR2 is non- vanishing and proportional to R, namely 5R2 = —Rg(t). 
However, we cannot find the full solution perturbatively. A further exact solution can also 
be found and reads 

R p (t,r) = ^l . (56) 

In this case, the equation for R(t) depends on the equation of state. Furthermore, imposing 
the condition that p r = we recover the scale factor of the flat universe, namely R(t) ~ ti, 
while requiring pg = we recover the result of the closed universe, namely R(t) ~ 9 — sin#. 

The structure of the small-time perturbative expansion in the open universe is such that 
new terms corresponding to new growing modes are introduced into the recursive relations. 
For the flat universe we clearly find a dominant term in t, of the type 5R^(t) ~ ^at+i 
for large values of N, foreseeing a logarithmic correction to the scale factor. Moreover, the 

2 

correction term itself grows at relatively small values of time. Since R(t) = At 3, the Hubble 
constant is proportional to A 3 / 2 , and we obtain 5Rn ~ A N+l t ^2 N ti , an d Sn ~ A N r \l, +1 . 



The higher-order terms are, therefore, important for t» ~ Ar 2 ~ Hor 3 , in the case of the 
flat universe. 



VII. CONCLUSIONS 

We analyzed the perturbation of a homogeneous Friedmann-Robertson- Walker space- 
time as a power series in the inverse of the radial parameter for the S-wave modes. The 
perturbative solutions depend not only on whether the universe is open, closed or flat, but 
also on the epoch when the perturbation starts, namely if it starts at early or at late times 
compared to the age of the universe. The standard expressions for the density contrast in 
both linear and non-linear regimes are contained in our results. 
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For the open and closed universes the full expression for the density contrast cannot be 
obtained, since one needs to solve a set of inhomogeneous recursive differential equations 
for the perturbation functions, where at a given order the lower-order solutions become the 
inhomogeneous part of the differential equation. However, the homogeneous part of the 
equations can be solved in closed form in all cases, and contains, in addition to the usual 
modes, infinitely many growing modes. The inhomogeneous part of the equation contains 
further valuable informations as far as the nonlinear modes are concerned. We have solved 
the full inhomogenous equation for the open matter-dominated universe in two limits: small 
and large times. At small times, the open universe behaves as a flat universe, with the 
difference that the nonlinear growing modes appear at lower orders in the perturbation 
series. The analysis of the large-time solutions, which covers both radiation and matter- 
dominated eras, indicates the disappearance of the growing modes at future times. 

For the flat universe, the differential equation simplifies since the scale factor depends 
explicitly on time and not parametrically through trigonometric expressions as in the closed 
and open universes. We have solved the inhomogeneous differential equation for the flat 
universe and have obtained the full expression for the linear and nonlinear growing modes 
(which fully agrees with the standard results). These modes grow faster as one goes to 
higher orders in the perturbation series. This points towards a bottom-up scenario for the 
structure formation and a universe dominated by cold dark matters. 

In order to verify the observational consequences of the modes growing faster than the 
standard modes of the linear regime, it is necessary to rewrite the results in terms of ob- 
servable quantities such as the luminosity distance rather than the variables (t,r). This is 
presently under investigation fl6fl . 

We thank F. Brandt for his help with the Maple programs. R.M. thanks A. Albrecht 
for helpful discussions. This work was supported by Conselho Federal de Desenvolvimento 
Cientifico e Tecnologico (CNPq-Brazil) and Fundagao de Amparo a Pesquisa do Estado de 
Sao Paulo (FAPESP). 
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